A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

A chromosome-level assembly of the seed beetle Callosobruchus maculatus genome with annotation of its repetitive elements.

Callosobruchus maculatus is a major agricultural pest of legume crops worldwide and an established model system in ecology and evolution. Yet, current molecular biological resources for this species are limited. Here, we employ Hi-C sequencing to generate a greatly improved genome assembly and we annotate its repetitive elements in a dedicated in-depth effort where we manually curate and classify the most abundant unclassified repeat subfamilies. We present a scaffolded chromosome-level assembly, which is 1.01 Gb in total length with 86% being contained within the 9 autosomes and the X chromosome. Repetitive sequences accounted for 70% of the total assembly. DNA transposons covered 18% of the genome, with the most abundant superfamily being Tc1-Mariner (9.75% of the genome). This new chromosome-level genome assembly of C. maculatus will enable future genetic and evolutionary studies not only of this important species but of beetles more generally.

B cell polygenic risk scores associate with anti-dsDNA antibodies and nephritis in systemic lupus erythematosus.

OBJECTIVE: B cell function and autoantibodies are important in SLE pathogenesis. In this work, we aimed to investigate the impact of cumulative SLE B cell genetics on SLE subphenotype and autoantibody profile. METHODS: Female patients with SLE (n=1248) and healthy controls (n=400) were genotyped using Illumina's Global Screening Array. Two polygenic risk scores (PRSs), one representing B cell genes and the other B cell activation genes, were calculated for each individual using risk loci for SLE in genes assigned to B cell-related pathways according to the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology and Reactome Databases. RESULTS: Double-stranded DNA (dsDNA) antibodies were more prevalent among patients with a high compared with a low SLE B cell PRS (OR 1.47 (1.07 to 2.01), p=0.018), and effect sizes were augmented in patients with human leucocyte antigen (HLA) risk haplotypes HLA-DRB1*03:01 and HLA-DRB1*15:01 (DRB1*03/15 -/- (OR 0.99 (0.56 to 1.77), p=0.98; DRB1*03/15 +/- or -/+ (OR 1.64 (1.06 to 2.54), p=0.028; and DRB1*03/15 +/+ (OR 4.47 (1.21 to 16.47), p=0.024). Further, a high compared with a low B cell PRS was associated with low complement levels in DRB1*03/15 +/+ patients (OR 3.92 (1.22 to 12.64), p=0.022). The prevalence of lupus nephritis (LN) was higher in patients with a B cell activation PRS above the third quartile compared with patients below (OR 1.32 (1.00 to 1.74), p=0.048). CONCLUSIONS: High genetic burden related to B cell function is associated with dsDNA antibody development and LN. Assessing B cell PRSs may be important in order to determine immunological pathways influencing SLE and to predict clinical phenotype.

Large-scale phage-based screening reveals extensive pan-viral mimicry of host short linear motifs.

Viruses mimic host short linear motifs (SLiMs) to hijack and deregulate cellular functions. Studies of motif-mediated interactions therefore provide insight into virus-host dependencies, and reveal targets for therapeutic intervention. Here, we describe the pan-viral discovery of 1712 SLiM-based virus-host interactions using a phage peptidome tiling the intrinsically disordered protein regions of 229 RNA viruses. We find mimicry of host SLiMs to be a ubiquitous viral strategy, reveal novel host proteins hijacked by viruses, and identify cellular pathways frequently deregulated by viral motif mimicry. Using structural and biophysical analyses, we show that viral mimicry-based interactions have similar binding strength and bound conformations as endogenous interactions. Finally, we establish polyadenylate-binding protein 1 as a potential target for broad-spectrum antiviral agent development. Our platform enables rapid discovery of mechanisms of viral interference and the identification of potential therapeutic targets which can aid in combating future epidemics and pandemics.

Experimental Life History Evolution Results in Sex-specific Evolution of Gene Expression in Seed Beetles.

The patterns of reproductive timing and senescence vary within and across species owing to differences in reproductive strategies, but our understanding of the molecular underpinnings of such variation is incomplete. This is perhaps particularly true for sex differences. We investigated the evolution of sex-specific gene expression associated with life history divergence in replicated populations of the seed beetle Acanthoscelides obtectus, experimentally evolving under (E)arly or (L)ate life reproduction for >200 generations which has resulted in strongly divergent life histories. We detected 1,646 genes that were differentially expressed in E and L lines, consistent with a highly polygenic basis of life history evolution. Only 30% of differentially expressed genes were similarly affected in males and females. The evolution of long life was associated with significantly reduced sex differences in expression, especially in non-reproductive tissues. The expression differences were overall more pronounced in females, in accordance with their greater phenotypic divergence in lifespan. Functional enrichment analysis revealed differences between E and L beetles in gene categories previously implicated in aging, such as mitochondrial function and defense response. The results show that divergent life history evolution can be associated with profound changes in gene expression that alter the transcriptome in a sex-specific way, highlighting the importance of understanding the mechanisms of aging in each sex.

A possible genomic footprint of polygenic adaptation on population divergence in seed beetles?

Efforts to unravel the genomic basis of incipient speciation are hampered by a mismatch between our toolkit and our understanding of the ecology and genetics of adaptation. While the former is focused on detecting selective sweeps involving few independently acting or linked speciation genes, the latter states that divergence typically occurs in polygenic traits under stabilizing selection. Here, we ask whether a role of stabilizing selection on polygenic traits in population divergence may be unveiled by using a phenotypically informed integrative approach, based on genome-wide variation segregating in divergent populations. We compare three divergent populations of seed beetles (Callosobruchus maculatus) where previous work has demonstrated a prominent role for stabilizing selection on, and population divergence in, key life history traits that reflect rate-dependent metabolic processes. We derive and assess predictions regarding the expected pattern of covariation between genetic variation segregating within populations and genetic differentiation between populations. Population differentiation was considerable (mean F ST = 0.23-0.26) and was primarily built by genes showing high selective constraints and an imbalance in inferred selection in different populations (positive Tajima's D NS in one and negative in one), and this set of genes was enriched with genes with a metabolic function. Repeatability of relative population differentiation was low at the level of individual genes but higher at the level of broad functional classes, again spotlighting metabolic genes. Absolute differentiation (d XY) showed a very different general pattern at this scale of divergence, more consistent with an important role for genetic drift. Although our exploration is consistent with stabilizing selection on polygenic metabolic phenotypes as an important engine of genome-wide relative population divergence and incipient speciation in our study system, we note that it is exceedingly difficult to firmly exclude other scenarios.

Proteome-scale mapping of binding sites in the unstructured regions of the human proteome.

Specific protein-protein interactions are central to all processes that underlie cell physiology. Numerous studies have together identified hundreds of thousands of human protein-protein interactions. However, many interactions remain to be discovered, and low affinity, conditional, and cell type-specific interactions are likely to be disproportionately underrepresented. Here, we describe an optimized proteomic peptide-phage display library that tiles all disordered regions of the human proteome and allows the screening of ~ 1,000,000 overlapping peptides in a single binding assay. We define guidelines for processing, filtering, and ranking the results and provide PepTools, a toolkit to annotate the identified hits. We uncovered >2,000 interaction pairs for 35 known short linear motif (SLiM)-binding domains and confirmed the quality of the produced data by complementary biophysical or cell-based assays. Finally, we show how the amino acid resolution-binding site information can be used to pinpoint functionally important disease mutations and phosphorylation events in intrinsically disordered regions of the proteome. The optimized human disorderome library paired with PepTools represents a powerful pipeline for unbiased proteome-wide discovery of SLiM-based interactions.

In silico promoter analysis and functional validation identify CmZFH, the co-regulator of hypoxia-responsive genes CmScylla and CmLPCAT.

Oxygen (O2) plays an essential role in aerobic organisms including terrestrial insects. Under hypoxic stress, the cowpea bruchid (Callosobruchus maculatus) ceases feeding and growth. However, larvae, particularly 4th instar larvae exhibit very high tolerance to hypoxia and can recover normal growth once brought to normoxia. To better understand the molecular mechanism that enables insects to cope with low O2 stress, we performed RNA-seq to distinguish hypoxia-responsive genes in midguts and subsequently identified potential common cis-elements in promoters of hypoxia-induced and -repressed genes, respectively. Selected elements were subjected to gel-shift and transient transfection assays to confirm their cis-regulatory function. Of these putative common cis-elements, AREB6 appeared to regulate the expression of CmLPCAT and CmScylla, two hypoxia-induced genes. CmZFH, the putative AREB6-binding protein, was hypoxia-inducible. Transient expression of CmZFH in Drosophila S2 cells activated CmLPCAT and CmScylla, and their induction was likely through interaction of CmZFH with AREB6. Binding to AREB6 was further confirmed by bacterially expressed CmZFH recombinant protein. Deletion analyses indicated that the N-terminal zinc-finger cluster of CmZFH was the key AREB6-binding domain. Through in silico and experimental exploration, we discovered novel transcriptional regulatory components associated with gene expression dynamics under hypoxia that facilitated insect survival.

Large scale discovery of coronavirus-host factor protein interaction motifs reveals SARS-CoV-2 specific mechanisms and vulnerabilities.

Viral proteins make extensive use of short peptide interaction motifs to hijack cellular host factors. However, most current large-scale methods do not identify this important class of protein-protein interactions. Uncovering peptide mediated interactions provides both a molecular understanding of viral interactions with their host and the foundation for developing novel antiviral reagents. Here we describe a viral peptide discovery approach covering 23 coronavirus strains that provides high resolution information on direct virus-host interactions. We identify 269 peptide-based interactions for 18 coronaviruses including a specific interaction between the human G3BP1/2 proteins and an ΦxFG peptide motif in the SARS-CoV-2 nucleocapsid (N) protein. This interaction supports viral replication and through its ΦxFG motif N rewires the G3BP1/2 interactome to disrupt stress granules. A peptide-based inhibitor disrupting the G3BP1/2-N interaction dampened SARS-CoV-2 infection showing that our results can be directly translated into novel specific antiviral reagents.

Sexual conflict drives micro- and macroevolution of sexual dimorphism in immunity.

BACKGROUND: Sexual dimorphism in immunity is believed to reflect sex differences in reproductive strategies and trade-offs between competing life history demands. Sexual selection can have major effects on mating rates and sex-specific costs of mating and may thereby influence sex differences in immunity as well as associated host-pathogen dynamics. Yet, experimental evidence linking the mating system to evolved sexual dimorphism in immunity are scarce and the direct effects of mating rate on immunity are not well established. Here, we use transcriptomic analyses, experimental evolution and phylogenetic comparative methods to study the association between the mating system and sexual dimorphism in immunity in seed beetles, where mating causes internal injuries in females. RESULTS: We demonstrate that female phenoloxidase (PO) activity, involved in wound healing and defence against parasitic infections, is elevated relative to males. This difference is accompanied by concomitant sex differences in the expression of genes in the prophenoloxidase activating cascade. We document substantial phenotypic plasticity in female PO activity in response to mating and show that experimental evolution under enforced monogamy (resulting in low remating rates and reduced sexual conflict relative to natural polygamy) rapidly decreases female (but not male) PO activity. Moreover, monogamous females had evolved increased tolerance to bacterial infection unrelated to mating, implying that female responses to costly mating may trade off with other aspects of immune defence, an hypothesis which broadly accords with the documented sex differences in gene expression. Finally, female (but not male) PO activity shows correlated evolution with the perceived harmfulness of male genitalia across 12 species of seed beetles, suggesting that sexual conflict has a significant influence on sexual dimorphisms in immunity in this group of insects. CONCLUSIONS: Our study provides insights into the links between sexual conflict and sexual dimorphism in immunity and suggests that selection pressures moulded by mating interactions can lead to a sex-specific mosaic of immune responses with important implications for host-pathogen dynamics in sexually reproducing organisms.

An experimental test of temperature-dependent selection on mitochondrial haplotypes in Callosobruchus maculatus seed beetles.

Mitochondrial DNA (mtDNA) consists of few but vital maternally inherited genes that interact closely with nuclear genes to produce cellular energy. How important mtDNA polymorphism is for adaptation is still unclear. The assumption in population genetic studies is often that segregating mtDNA variation is selectively neutral. This contrasts with empirical observations of mtDNA haplotypes affecting fitness-related traits and thermal sensitivity, and latitudinal clines in mtDNA haplotype frequencies. Here, we experimentally test whether ambient temperature affects selection on mtDNA variation, and whether such thermal effects are influenced by intergenomic epistasis due to interactions between mitochondrial and nuclear genes, using replicated experimental evolution in Callosobruchus maculatus seed beetle populations seeded with a mixture of different mtDNA haplotypes. We also test for sex-specific consequences of mtDNA evolution on reproductive success, given that mtDNA mutations can have sexually antagonistic fitness effects. Our results demonstrate natural selection on mtDNA haplotypes, with some support for thermal environment influencing mtDNA evolution through mitonuclear epistasis. The changes in male and female reproductive fitness were both aligned with changes in mtDNA haplotype frequencies, suggesting that natural selection on mtDNA is sexually concordant in stressful thermal environments. We discuss the implications of our findings for the evolution of mtDNA.

Negative frequency dependent selection contributes to the maintenance of a global polymorphism in mitochondrial DNA.

BACKGROUND: Understanding the forces that maintain diversity across a range of scales is at the very heart of biology. Frequency-dependent processes are generally recognized as the most central process for the maintenance of ecological diversity. The same is, however, not generally true for genetic diversity. Negative frequency dependent selection, where rare genotypes have an advantage, is often regarded as a relatively weak force in maintaining genetic variation in life history traits because recombination disassociates alleles across many genes. Yet, many regions of the genome show low rates of recombination and genetic variation in such regions (i.e., supergenes) may in theory be upheld by frequency dependent selection. RESULTS: We studied what is essentially a ubiquitous life history supergene (i.e., mitochondrial DNA) in the fruit fly Drosophila subobscura, showing sympatric polymorphism with two main mtDNA genotypes co-occurring in populations world-wide. Using an experimental evolution approach involving manipulations of genotype starting frequencies, we show that negative frequency dependent selection indeed acts to maintain genetic variation in this region. Moreover, the strength of selection was affected by food resource conditions. CONCLUSIONS: Our work provides novel experimental support for the view that balancing selection through negative frequency dependency acts to maintain genetic variation in life history genes. We suggest that the emergence of negative frequency dependent selection on mtDNA is symptomatic of the fundamental link between ecological processes related to resource use and the maintenance of genetic variation.

The genomic footprint of sexual conflict.

Genes with sex-biased expression show a number of unique properties and this has been seen as evidence for conflicting selection pressures in males and females, forming a genetic 'tug-of-war' between the sexes. However, we lack studies of taxa where an understanding of conflicting phenotypic selection in the sexes has been linked with studies of genomic signatures of sexual conflict. Here, we provide such a link. We used an insect where sexual conflict is unusually well understood, the seed beetle Callosobruchus maculatus, to test for molecular genetic signals of sexual conflict across genes with varying degrees of sex-bias in expression. We sequenced, assembled and annotated its genome and performed population resequencing of three divergent populations. Sex-biased genes showed increased levels of genetic diversity and bore a remarkably clear footprint of relaxed purifying selection. Yet, segregating genetic variation was also affected by balancing selection in weakly female-biased genes, while male-biased genes showed signs of overall purifying selection. Female-biased genes contributed disproportionally to shared polymorphism across populations, while male-biased genes, male seminal fluid protein genes and sex-linked genes did not. Genes showing genomic signatures consistent with sexual conflict generally matched life-history phenotypes known to experience sexually antagonistic selection in this species. Our results highlight metabolic and reproductive processes, confirming the key role of general life-history traits in sexual conflict.

Identification of novel ejaculate proteins in a seed beetle and division of labour across male accessory reproductive glands.

The male ejaculate contains a multitude of seminal fluid proteins (SFPs), many of which are key reproductive molecules, as well as sperm. However, the identification of SFPs is notoriously difficult and a detailed understanding of this complex phenotype has only been achieved in a few model species. We employed a recently developed proteomic method involving whole-organism stable isotope labelling coupled with proteomic and transcriptomic analyses to characterize ejaculate proteins in the seed beetle Callosobruchus maculatus. We identified 317 proteins that were transferred to females at mating, and a great majority of these showed signals of secretion and were highly male-biased in expression in the abdomen. These male-derived proteins were enriched with proteins involved in general metabolic and catabolic processes but also with proteolytic enzymes and proteins involved in protection against oxidative stress. Thirty-seven proteins showed significant homology with SFPs previously identified in other insects. However, no less than 92 C. maculatus ejaculate proteins were entirely novel, receiving no significant blast hits and lacking homologs in extant data bases, consistent with a rapid and divergent evolution of SFPs. We used 3D micro-tomography in conjunction with proteomic methods to identify 5 distinct pairs of male accessory reproductive glands and to show that certain ejaculate proteins were only recovered in certain male glands. Finally, we provide a tentative list of 231 candidate female-derived reproductive proteins, some of which are likely important in ejaculate processing and/or sperm storage.

The Evolution of Dark Matter in the Mitogenome of Seed Beetles.

Animal mitogenomes are generally thought of as being economic and optimized for rapid replication and transcription. We use long-read sequencing technology to assemble the remarkable mitogenomes of four species of seed beetles. These are the largest circular mitogenomes ever assembled in insects, ranging from 24,496 to 26,613 bp in total length, and are exceptional in that some 40% consists of non-coding DNA. The size expansion is due to two very long intergenic spacers (LIGSs), rich in tandem repeats. The two LIGSs are present in all species but vary greatly in length (114-10,408 bp), show very low sequence similarity, divergent tandem repeat motifs, a very high AT content and concerted length evolution. The LIGSs have been retained for at least some 45 my but must have undergone repeated reductions and expansions, despite strong purifying selection on protein coding mtDNA genes. The LIGSs are located in two intergenic sites where a few recent studies of insects have also reported shorter LIGSs (>200 bp). These sites may represent spaces that tolerate neutral repeat array expansions or, alternatively, the LIGSs may function to allow a more economic translational machinery. Mitochondrial respiration in adult seed beetles is based almost exclusively on fatty acids, which reduces the need for building complex I of the oxidative phosphorylation pathway (NADH dehydrogenase). One possibility is thus that the LIGSs may allow depressed transcription of NAD genes. RNA sequencing showed that LIGSs are partly transcribed and transcriptional profiling suggested that all seven mtDNA NAD genes indeed show low levels of transcription and co-regulation of transcription across sexes and tissues.

Mating Changes Sexually Dimorphic Gene Expression in the Seed Beetle Callosobruchus maculatus.

Sexually dimorphic phenotypes arise largely from sex-specific gene expression, which has mainly been characterized in sexually naïve adults. However, we expect sexual dimorphism in transcription to be dynamic and dependent on factors such as reproductive status. Mating induces many behavioral and physiological changes distinct to each sex and is therefore expected to activate regulatory changes in many sex-biased genes. Here, we first characterized sexual dimorphism in gene expression in Callosobruchus maculatus seed beetles. We then examined how females and males respond to mating and how it affects sex-biased expression, both in sex-limited (abdomen) and sex-shared (head and thorax) tissues. Mating responses were largely sex-specific and, as expected, females showed more genes responding compared with males (∼2,000 vs. ∼300 genes in the abdomen, ∼500 vs. ∼400 in the head and thorax, respectively). Of the sex-biased genes present in virgins, 16% (1,041 genes) in the abdomen and 17% (243 genes) in the head and thorax altered their relative expression between the sexes as a result of mating. Sex-bias status changed in 2% of the genes in the abdomen and 4% in the head and thorax following mating. Mating responses involved de-feminization of females and, to a lesser extent, de-masculinization of males relative to their virgin state: mating decreased rather than increased dimorphic expression of sex-biased genes. The fact that regulatory changes of both types of sex-biased genes occurred in both sexes suggests that male- and female-specific selection is not restricted to male- and female-biased genes, respectively, as is sometimes assumed.

Divergent evolution of life span associated with mitochondrial DNA evolution.

Mitochondria play a key role in ageing. The pursuit of genes that regulate variation in life span and ageing have shown that several nuclear-encoded mitochondrial genes are important. However, the role of mitochondrial encoded genes (mtDNA) is more controversial and our appreciation of the role of mtDNA for the evolution of life span is limited. We use replicated lines of seed beetles that have been artificially selected for long or short life for >190 generations, now showing dramatic phenotypic differences, to test for a possible role of mtDNA in the divergent evolution of ageing and life span. We show that these divergent selection regimes led to the evolution of significantly different mtDNA haplotype frequencies. Selection for a long life and late reproduction generated positive selection for one specific haplotype, which was fixed in most such lines. In contrast, selection for reproduction early in life led to both positive selection as well as negative frequency-dependent selection on two different haplotypes, which were both present in all such lines. Our findings suggest that the evolution of life span was in part mediated by mtDNA, providing support for the emerging general tenet that adaptive evolution of life-history syndromes may involve mtDNA.

The De Novo Transcriptome and Its Functional Annotation in the Seed Beetle Callosobruchus maculatus.

Despite their unparalleled biodiversity, the genomic resources available for beetles (Coleoptera) remain relatively scarce. We present an integrative and high quality annotated transcriptome of the beetle Callosobruchus maculatus, an important and cosmopolitan agricultural pest as well as an emerging model species in ecology and evolutionary biology. Using Illumina sequencing technology, we sequenced 492 million read pairs generated from 51 samples of different developmental stages (larvae, pupae and adults) of C. maculatus. Reads were de novo assembled using the Trinity software, into a single combined assembly as well as into three separate assemblies based on data from the different developmental stages. The combined assembly generated 218,192 transcripts and 145,883 putative genes. Putative genes were annotated with the Blast2GO software and the Trinotate pipeline. In total, 33,216 putative genes were successfully annotated using Blastx against the Nr (non-redundant) database and 13,382 were assigned to 34,100 Gene Ontology (GO) terms. We classified 5,475 putative genes into Clusters of Orthologous Groups (COG) and 116 metabolic pathways maps were predicted based on the annotation. Our analyses suggested that the transcriptional specificity increases with ontogeny. For example, out of 33,216 annotated putative genes, 51 were only expressed in larvae, 63 only in pupae and 171 only in adults. Our study illustrates the importance of including samples from several developmental stages when the aim is to provide an integrative and high quality annotated transcriptome. Our results will represent an invaluable resource for those working with the ecology, evolution and pest control of C. maculatus, as well for comparative studies of the transcriptomics and genomics of beetles more generally.

Negative frequency dependent selection on sympatric mtDNA haplotypes in Drosophila subobscura.

BACKGROUND: Recent experimental evidence for selection on mitochondrial DNA (mtDNA) has prompted the question as to what processes act to maintain within-population variation in mtDNA. Balancing selection though negative frequency dependent selection (NFDS) among sympatric haplotypes is a possibility, but direct empirical evidence for this is very scarce. FINDINGS: We extend the previous findings of a multi-generation replicated cage experiment in Drosophila subobscura, where mtDNA polymorphism was maintained in a laboratory setting. First, we use a set of Monte Carlo simulations to show that the haplotype frequency dynamics observed are inconsistent with genetic drift alone and most closely match those expected under NFDS. Second, we show that haplotype frequency changes over time were significantly different from those expected under either genetic drift or positive selection but were consistent with those expected under NFSD. CONCLUSIONS: Collectively, our analyses provide novel support for NFDS on mtDNA haplotypes, suggesting that mtDNA polymorphism may at least in part be maintained by balancing selection also in natural populations. We very briefly discuss the possible mechanisms that might be involved.

Genome size correlates with reproductive fitness in seed beetles.

The ultimate cause of genome size (GS) evolution in eukaryotes remains a major and unresolved puzzle in evolutionary biology. Large-scale comparative studies have failed to find consistent correlations between GS and organismal properties, resulting in the 'C-value paradox'. Current hypotheses for the evolution of GS are based either on the balance between mutational events and drift or on natural selection acting upon standing genetic variation in GS. It is, however, currently very difficult to evaluate the role of selection because within-species studies that relate variation in life-history traits to variation in GS are very rare. Here, we report phylogenetic comparative analyses of GS evolution in seed beetles at two distinct taxonomic scales, which combines replicated estimation of GS with experimental assays of life-history traits and reproductive fitness. GS showed rapid and bidirectional evolution across species, but did not show correlated evolution with any of several indices of the relative importance of genetic drift. Within a single species, GS varied by 4-5% across populations and showed positive correlated evolution with independent estimates of male and female reproductive fitness. Collectively, the phylogenetic pattern of GS diversification across and within species in conjunction with the pattern of correlated evolution between GS and fitness provide novel support for the tenet that natural selection plays a key role in shaping GS evolution.